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Abstract 



The sensitivity of nucleon-nucleus elastic scattering to the off-shell behavior 
of realistic nucleon-nucleon interactions is investigated when on-shell equiva- 
lent nucleon-nucleon potentials are used. The study is based on applications 
of the full-folding optical model potential for an explicit treatment of the off- 
shell behavior of the nucleon-nucleon effective interaction. Applications were 
made at beam energies between 40 and 500 MeV for proton scattering from 
40 Ca and 208 Pb. We use the momentum-dependent Paris potential and its lo- 
cal on-shell equivalent as obtained with the Gelfand-Levitan and Marchenko 
inversion formalism for the two nucleon Schrodinger equation. Full-folding 
calculations for nucleon-nucleus scattering show small fluctuations in the cor- 
responding observables. This implies that off-shell features of the NN inter- 
action cannot be unambiguously identified with these processes. Inversion 
potentials were also constructed directly from NN phase-shift data (SM94) 
in the 0-1.3 GeV energy range. Their use in proton-nucleus scattering above 
200 MeV provide a superior description of the observables relative to those 
obtained from current realistic NN potentials. Limitations and scope of our 
findings are presented and discussed. 

PACS numbers: 24.10.-i, 24.10.Ht, 13.75.Cs, 24.80. +y, 25.40. Cm 
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I. INTRODUCTION 



Theoretical studies of the two-nucleon interaction in their off-shell domain have a long 
standing tradition |IJ and this topic is of renewed interest for designated experiments at sev- 
eral accelerators laboratories. Few- and many-body systems offer generally the possibilities 
for such studies but not seldom have such endeavors ended prematurely due to lack of statis- 
tics in data taking or an incomplete and inconclusive theory. Current interest on this issue 
comes from theory groups who have independently developed NN potentials which account 
reasonable well for the two-body phase-shifts at energies below pion production threshold 
These potentials are manifestly different in aspects such as their off-shell behavior 
and the device of an experimental discrimination among them would be very valuable for a 
comprehensive understanding of particularities for each of the currently acceptable models. 

In the past much hope was given to NN bremsstrahlung since this three body reaction 
is within a Born approximation theoretically well defined As a result of this simple 

reaction mechanism one obtains a good perspective of the link between the (NNj) observ- 
ables and half-off-shell NN t matrices. Using the available theoretical developments in the 
field of bremsstrahlung it is now known how to explain the available data with different NN 
potentials as far as they reproduce the NN on-shell data very well. This result has been 
independently validated and remains a surprise in view of the obvious differences in the 
half-off-shell t matrices from different model potentials. In other words, the expectation 
that on-shell equivalent NN potentials would provide distinctive bremsstrahlung predictions 
was disillusioned with more complete and reliable calculations. 

It can be argued that bremsstrahlung involves only half-off-shell t matrices and thus 
this reaction is rather confined in phase-space. Microscopic models of nucleon-nucleus 
(NA) scattering in a full-folding framework ]7|J| do not suffer from such limitations as the 
optical potential depends explicitly upon the effective interaction fully off-shell and therefore 
constitute a wider frame for investigating the NN interaction off-shell. 

Recently, significant advances have been made in accurately handling the off-shell degrees 
of freedom in NA elastic scattering 0-12|. Here, studies have demonstrated that an accurate 



treatment of the off-shell behavior of the NN interaction is needed for a proper account of 
the theory. Irrespective of obvious improvements in describing the data at projectile energies 
below 400 MeV, the calculations still show some systematic deficiencies. In particular the 
misfit of spin observables at small momentum transfers, < lfm -1 , is not understood. The 
origin of such discrepancies with the data could be attributed to intrinsic limitations of the 
bare NN potential with respect to the phenomenology (especially at the higher energies), to 
the simplifications implicit in the model for the NN effective interaction, or to the fact that 
the optical model potential has only been developed to its lowest order. Above 400 MeV, 
full folding model results quickly deteriorate as they fail to properly describe the observables 
0- 

Despite the significant advances in full-folding model calculations it remains difficult 
to identify, at the level of the scattering observables, distinctive off-shell features of NN 
interactions. The difficulty is mainly attributed to non negligible differences on-shell among 
currently available realistic NN potentials [ |i~3H . 

In this paper we investigate the sensitivity of NA scattering when using bare NN poten- 
tials which are equivalent on-shell. For that purpose we have chosen the Paris potential (i.) 



2 



in its full momentum dependent form and (ii.) a local equivalent potential generated with 
Gelfand-Levitan-Marchenko inversion from the Paris potential phase-shifts and deuteron 
bound state. The inversion potential yields different off-shell continuations as the origi- 
nal Paris potential. These differences and associated correlations are investigated in the 
full-folding NA optical model. 

We find for these two alternatives differences of about 5% overall in NA observables and 
this rises doubts that this sensitivity could be used to infer off-shell peculiarities of NN 
interactions in general. Based on this result we conclude that quantum inversion provides 
a practical and accurate connection between NN and NA elastic scattering. As mediating 
formalism may Gelfand-Levitan-Marchenko quantum inversion use only experimental NN 
data as input and thus link model-independent two-body with many-body data. 

This article is organized as follows. In section II we outline the theoretical background for 
the present work. We introduce and discuss aspects and assumptions implicit in the Gelfand- 
Levitan and Marchenko inversion method as well as the in-medium full-folding model of 
the optical potential for NA elastic scattering. In section III we discuss the sensitivity of 
NA scattering by calculating the scattering observables from the Paris potential and from 
an inversion NN potential constructed from the Paris NN phase-shifts. In this way, our 
analysis becomes explicitly dependent only from the off-shell differences between the two 
potentials. In section IV we construct inversion NN potentials directly from the NN data, 
including approximately the NN phase-shifts above pion production threshold as a way to 
guide the intermediate energy properties of the potential. Calculations for proton-nucleus 
elastic scattering in the 40-500 MeV range and from different targets are also discussed. 
Finally, in section V we present a summary and draw conclusions from our work. 

II. THEORETICAL BACKGROUND 

A. Two Nucleon Potentials from quantum inversion 

The starting point for the NN interaction is generally a relativistic scattering description 
and it is beyond any doubt, that a realistic potential must show strong angular momentum 
and spin dependencies and thus channel dependence. An explicit momentum dependence 
or nonlocality is predicted in all relativistic potential models. Allowance for the relativistic 
nature of the problem in these models is in the best cases limited to relativistic kinematics 
and simplified calculation of selected exchange diagrams. The full dynamics of potentials 
requires a derivation using quantum field theory. This not only ensures Lorentz covariance 
but shows the way how the actual potential must be constructed. For this purpose, different 
formalism have been proposed and they determine the off-energy shell or off-mass shell am- 
plitudes. Usually the derivation of such relativistic equations also determines the propaga- 
tors and the potential which contains all information about the dynamics of the interaction. 
An accepted four-dimensional formulation of quantum field theory is the covariant Bethe- 
Salpeter equation and the three-dimensional relativistic equations are obtained on the basis 
of quasi-potential methods. Another field theoretical methods use a generalization of the 
Schrodinger equation in the form of the Tomonaga-Schwinger equation and in the form of 
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a covariant Hamiltonian formulation. Irrespective of their foundation and motivation, the 
field theoretical models are using more or less phenomenology to fit data. 

Well fitted and widely applied within these models are the boson exchange potentials. 



Representative within the boson exchange models are the Nijmegen |L4j and Bonn potentials 
0. The Paris potential also belongs to this category and describes the long- and medium- 
range interaction with single and correlated two-pion exchange and heavier meson exchanges 
0. It fits the experimental two nucleon data reasonably well and has tradition in calculations 
of microscopic optical model potentials for NA scattering at low and medium energies. It is 
generally used in its Yukawa parameterization including an explicit momentum dependent 
term @ 

V = V a + -p 2 V b + V b -p 2 , (1) 
m m 

where m is the nucleon mass and p the relative momentum operator. This momentum 
dependence emulates a hard core repulsion and produces a divergence of the phase-shift 
5(k) as k — > oo. When compared with data, the Paris potential fit quite well the phase shift 
at energies up to 280 MeV. Above this energy we notice a rapidly increasing divergence. 
This characteristic strong repulsion at short distances is easy to identify in cross section and 



spin observables of NA scattering and reactions [15]. Despite this weakness we have chosen 
the Paris potential since we have confidence in the numerics of our Lippmann-Schwinger 
and Bethe-Goldstone calculations with it. 

With this investigation we aim to disclose effects which are caused by the momentum 
dependence and high energy phase shift discrepancies of the genuine Paris potential in off- 
shell t— and g— matrices. To this purpose we use the recently developed Gelfand-Levitan- 
Marchenko inversion of partial wave radial Schrodinger equations to generate phase equiva- 
lent local potentials to Paris potential phase shifts. These inversion algorithms distinguish 
inversion for single and coupled channels cases with and without a Coulomb reference poten- 



tial ||16|| . In other words, we obtain separately the hadronic part of the NN interaction from 
any set of np and pp phase shifts. To investigate the effect of different on-shell behaviors we 
also generate potentials from the latest experimental phase shift analysis SM94 by Arndt 
and collaborators Hl7j . 

The strong interaction inversion potentials V(r; LSJ,T; np) and V(r; LSJ,T; pp) are 
numerical solutions of Gelfand-Levitan or Marchenko integral equations, 

K(r, r') + F(r, r') + J K(r, s)F(s, r') ds = (2) 

and 

V(r) = ±-^-K(r,r), (3) 

for any specified radius and they are determined channel by channel. The input kernel 
F(x,y) is computed with the spectral information, Jost functions or S-matrices including 
deuteron binding energy and normalization constants. We use quantum inversion as trans- 
formation of given real phase shifts, which are specified within a finite energy interval 
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5(k) ~ 5(E) = {5(k)\E = [0,E max }\k = [0,k max ]}. (4) 

Thereafter they are smoothly extrapolated with a rational function which decays asymptot- 
ically @ 



5(E) ~ (J(A:) = (5(fc)|fc > fc max | lim ~ A;" 1 ) . (5) 

This transformation is unique for the class of potentials in which we are interested and which 
we assume physically significant. The resulting inversion potentials are real and local. 

It is evident that this inversion procedure produces a restricted phase equivalent potential 
to the Paris potential. The limited data input 5(E) for < E < E max is used to control 
the range of equivalence and the extrapolation thereafter to control the softness of the short 
range core interaction. 



B. In-medium full folding optical potential 

The optical potential for NA elastic scattering can be casted as the convolution of an anti- 
symmetrized effective interaction with the target ground state single particle wave functions 
PH|-P^,P[. In momentum space this one-body operator reads 

U(k',k;E) = J dp dp' ]T <Pi(p'){k',p' \F(E + e a )\k,p) A <f> a (P) (6) 

a<ep 

where E represents the energy of the incoming projectile and {4> a , e a } are the target ground 
state single-particle wave functions and corresponding energies. The momenta k(k') and 
pip') correspond to the initial (final) momenta of the projectile and target struck nucleon 
respectively. The two-nucleon interaction in Eq. accounts for multiple scattering of 
nucleons to all orders in the ladder approximation |T9| . |23|j . Although a general expression 
for this matrix can formally be defined, its practical implementation requires the device of 
a dynamical model for the effective two-nucleon interaction in the nucleus. The procedure 
we follow is that introduced in Ref. [0, where translational invariance properties of two- 
nucleon scattering in free space or infinite nuclear matter suggest the following ansatz for 
the two-body matrix 



k',p' | F(u) | k,p) = j^j dRe^<^ 



k' 



J Q+Q v ' ' 
— 1 — 



k) . (7) 



(27T) 

Here we have defined the initial and final two-nucleon center-of-mass (cm.) momenta, 

Q = k + p, Q' = k'+p', (8) 
and the corresponding relative momenta by 

K=l(k-p), K' = \(k'-p'). (9) 

The function {K,'\f^(uj; R)\k) corresponds to the matrix elements of a reduced two-body 
effective interaction. In the case of no dependence of the / matrix upon the spatial coordinate 
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(R) one restores total momentum conservation of the interacting pair as the radial integral 
in Eq. (^) leads to a cm. momentum conserving Dirac 5— function. 

A calculable expression for the optical potential in medium emerges after a systematic 
reduction of the many body propagator when represented in terms of the target ground state 
spectral function || . To lowest order in a series expansion of the two-body propagator in a 
finite nucleus this interaction can be identified with the g matrix solution of the Brueckner- 
Bethe-Goldstone equation for interacting nucleons in infinite nuclear matter evaluated at 
nuclear density p(R) in the nucleus. Furthermore, it becomes convenient at this point to 
substitute the single particle energies e a by an average value, e, and to use the Slater or 
Campi-Bouyssy approximations P||,PBJ to represent the ground state mixed density, i.e. 



p(p',p)= E « p /^^fl) /dp e to - p] } , (io) 

where p(R) is the local nuclear density at coordinate R and P represents the struck nucleon 
mean momentum defined by 

p = W+p')- (ii) 

The local momentum function k(R) sets the range of variation of the struck nucleon mean 
momenta upon collisions with the projectile and is obtained from the Slater or Campi- 
Bouyssy prescriptions. The local density function p(R) is defined in terms of the local 
momentum function by 

p{R) = 7^- 2 k 3 (R) • (12) 

With the above considerations the optical potential can be expressed in terms of the nuclear 
density and a Fermi averaged effective interaction obtained from interacting nuclear matter. 
This interaction retains nuclear medium correlations associated with the nuclear mean fields 
and Pauli blocking. The in-medium full-folding optical potential then reads 



(2vr)3 7 ^ 'p(R) 

x e[k(R)-P] {\{K - P - q)\g R+P {E + e;R)\l(K - P + q)) a . (13) 

Thus, the optical potential requires the calculation of g matrices off-shell as their relative 
momenta obey no constraints apart from those imposed by the ground state mixed density 
of the target. Furthermore, no assumptions are introduced on the nature of the momentum 
dependence of the optical potential, thus retaining all non localities arising from the genuine 
momentum dependence of the NN effective interaction and as prescribed by the full-folding 
integral. Actual calculations involve determining g matrices at several densities and over 
a wide range of total center-of-mass momenta, features fully accounted for in the present 
work. 

In the context of a medium independent internucleon interaction, as when the free t ma- 
trix is used to represent the NN effective interaction, the integral over the spatial coordinate 
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in Eq. (|T3|) can be performed separately from the motion of the target nucleons. With the 
use of Eq. (|10|) for the mixed density one recovers the expression for the full-folding optical 
potential in the zero density approximation J7j, namely U(k', k; E) — ► U (k', k; E), where 

U (k',k;E) = JdPp(P+lq,P-^{l(K-P-q)\t. +p (E + e)\l(K-P + q)) A . (14) 

The dependence of the optical potential on off-shell t— matrices becomes explicit in 
the above expression. The feasibility of the full-folding model to investigate particular 
signatures of the effective interaction off-shell will depend on the sensitivity of NA scattering 
observables to the use of t matrices with manifestly distinctive behaviors off-shell. An 
important constraint for such study is that the effective interactions, t matrices in this limit, 
agree on-shell. To the extent this constraint is met, one can attain the differences in the 
NA scattering observables to the differences of the interactions off-shell. 

A simple kinematical effect usually overlooked, but explicitly accounted for in our cal- 
culations, is that the full-folding approach calls for matrix elements of energy E + e in the 
laboratory frame. In the limit of the free t matrix for the NN interactions, the energy of the 
interacting pair in its cm. is given by E + e — 4-{K-\- P) 2 . Therefore, the maximum energy 
of the pair in its cm. is E + e, the energy of the beam plus the average binding energy of the 
target nucleons. This is to say that in the case of optical potentials for nucleons at 500 MeV, 
t matrices of up to ~1 GeV in the laboratory frame are required. This more demanding 
sampling of the NN effective interaction is a result of the unconstrained kinematics allowed 
by the Fermi motion of the nucleons in the nucleus. 

A few comments are pertinent regarding further approximations in the treatment of the 
t matrix which limit a clear assessment of the off-shell behavior of the NN interaction in 
NA scattering. A simplifying assumption, commonly used in some alternative full-folding 



calculations 10,|TT|| , is that the t matrix varies very weakly with respect to the NN cm. 
momentum K + P. Thus, the magnitude of this momentum is fixed to the (asymptotic) 
on-shell value of the incoming projectile, K Q , and the t matrix is approximated by 

( 15 ) 

Thus, t matrices are evaluated at a fixed energy equal in the NN cm. to one-half the 
energy of the beam. Here one has neglected all effects associated with the Fermi motion in 
the NN cm. momentum dependence. The resulting full-folding calculations samples the t 
matrix off-shell through its dependence on P in the relative momenta exclusively (see Eq. 
([14])). This approximation seems adequate at beam energies near 300 MeV. Its application at 
lower or higher energies, however, needs further considerations as the NN t matrix exhibits 
a sizable NN cm. momentum dependence ||. In the low energy region, apart from the 
fact that medium effects need to be incorporated in the model, the underlying kinematics 
prescribed by the full-folding yields the sampling of the t matrix in regions where it varies 
significantly as the low energy behavior of the interactions becomes dominant. In the high 
energy regime, in turn, difficulties arise from the opening of inelastic channels such as those 
associated with pion production or A— resonances. The actual merit of the theory of the 
optical potential needs to be assessed with a consistent incorporation of such additional 
degrees of freedom. 
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C. The NN effective interaction 



In the present approach, correlations associated with interacting nucleons in the nuclear 
medium are obtained from the NN effective interaction defined by the Brueckner-Bethe- 
Goldstone equation for symmetric nuclear matter. In momentum representation, the g 
matrix associated with interacting nucleons of total cm. momentum Q, starting energy u 
and nuclear density p satisfies 



(FZ'\V\k) + I dn" {k'\V\k") \l (Fl";uj;k F )(K" <?-(o;;P) 



K 



(16) 



with k F the nuclear matter Fermi momentum determined from the nuclear density via 



3tt 5 



1/3 



p{R) 



(17) 



The two-body propagator A models both Pauli blocking and the nuclear mean field effects 
in the propagation of intermediate states, 



A (q ; u ; k F ) 



Q(P + ;P_;k F ) 



uj + ii] - e(P + ; k F ) - e(P_; k F ) ' 



(18) 



with P± = ± q and Q the Pauli blocking function 



Q{P + ;P_;k F ) = e[ e(P + ; k F )-e F ]G[ e(P_; k F ) - e F ] . 



(19) 



Here the single particle energies e are defined in terms of the self-consistent nuclear matter 
fields, U vM) 



k 2 

e(k a ; k F ) = + Re [U NM (k a ; k F ) 



(20) 



where the mean fields U NM (k; k F ) are calculated self-consistently for the underlying bare 
NN interaction from 



U NM (k; k F )= (l(k-ka 



2 (/c k a )^ 



(21) 



Actual calculations of these mean fields have been made using the continuous prescription 
at the Fermi energy and simplifying the Pauli blocking function Q by its angle-averaged 
form. 



III. PARIS VERSUS PARIS INVERSION 

We have used the Paris potential to generate a set of NN phase-shifts which are taken 
as input to calculate the corresponding inversion potential. Thus, we make sure that Paris 
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and the inversion potential are equivalent on-shell within the accuracy obtained with the 
numeric algorithms used for implementing the quantum inversion method. The range of 
energies where the explicit phase-shifts are considered determine importantly the off-shell 
behavior of the inversion potential. This effect is observed when calculating, for example, 
the optical potential and observables for NA scattering from the free t matrix obtained from 
inversion potentials constructed from different sets of phase-shifts. Therefore, to ensure 
that properties of the inversion potential depend solely on the dynamical equations and 
not on the range of energies considered for the phase-shifts, we have verified that a set 
of Paris phase-shifts in the 0-1.3 GeV energy range is sufficient to construct an inversion 
potential that can be used up to 500 MeV in effective interaction calculations (Sec. II. C). 
Above 1.3 GeV, the inversion algorithm assumes a smooth decrease to zero at infinity for the 
phase-shifts (Eq. |5]). This extrapolation departs strongly from the Paris potential behavior. 

In Fig. [I] we show the Paris phase-shifts (crosses) and the phase-shifts obtained from 
its corresponding inversion potential (full line) for selected NN channels (L < 2) and in 
the 0-1.2 GeV energy range. In general, for most of the NN channels we have considered 
the agreement is excellent. Some differences are observed in the (coupled) 3 Di state above 
400 MeV. Altogether, we can conclude that the calculated inversion potential is phase- 
equivalent to the Paris potential. Their differences in a many-body system should come 
from the intrinsic properties of the two potentials and provided that the off-shell sampling 
is compatible with the energy range used to set the phase equivalence. 

The inversion potential differs on its off-shell content from the original Paris potential as 
the former is static and local, whereas the latter is momentum dependent. These differences 
and associated correlations to all orders are now investigated in the context of the full-folding 
model of the optical potential for NA scattering and as described in Sec. II. B. 



A. Sensitivity to off* shell effects in NA scattering 

We have calculated both t matrices and in-medium g matrices from the Paris and its 
inversion potentials. The inversion scheme was used for all the NN channels with total 
angular momentum J < 2. For channels with J > 2 the genuine Paris potential is used. The 
nucleon-nucleus optical potential was calculated using the g matrix as effective interaction 
since medium effects have been proved important even for nucleons with incident energies of 
400 MeV |J. The corresponding g matrices were calculated solving Eq. (16) using standard 



matrix inversion methods |27J . The R— dependence in the full-folding integral was obtained 
by calculating g matrices at different densities as obtained from different values for hp up 
to 1.4 fm . For an accurate off-shell sampling of the NN effective interaction, g matrices 
were calculated at several values of the total NN cm. momentum in the 0-7 fm -1 interval, 
with higher density of points in the region it varies most rapidly as a function of the cm. 
momentum. Contributions associated with the deuteron bound state singularity were also 
included ||. The ground state nuclear densities and average binding energies (e in Eq. fll3]) ) 
are the same as in Ref. 0. 

Calculations of differential cross sections (da/dfl) and analyzing powers (A y ) were made 
for proton elastic scattering from 40 Ca and 208 Pb in the 40-400 MeV energy range. Results 
for the spin rotation function have been omitted for brevity as they exhibit similar behavior 
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as those observed in A y . In Figs. ^| and |^ we present the calculated scattering observables, 
as function of the center of mass scattering angle (0 c . m .) in the 40-MeV application, and as 
a function of the momentum transfer (q) at 400 MeV. The data for p+ 40 Ca scattering data 
at 40 MeV were taken from Ref. [^,^]. In the case of p+ 208 Pb, the data at 400 MeV were 
taken from Ref. |30| . The full and dashed curves represent results for the Paris potential and 
corresponding inversion respectively. Results at different energies and for different targets 
show comparable quantitative differences. Also, very similar differences in the observables 
are observed when using a free t matrix as input in the full-folding calculations. 

Some conclusions can be drawn. The two phase-equivalent potentials give a qualitatively 
similar description of the scattering observables, regardless their very different intrinsic struc- 
ture. The differences between the two curves in Figs. ^| and || reflect the level of sensitivity 
to the off-shell behavior associated with each underlying NN effective interaction. Quanti- 
tatively, we observe up to a 10% difference in the magnitude of the calculated observables 
due to differences off-shell in the effective interaction, feature that is maintained in the en- 
ergy range under study. This weak sensitivity to differences off-shell of completely different 
NN bare potentials suggests that the full-folding model cannot unambiguously discriminate 
among on-shell equivalent NN potentials. Furthermore, this equivalence indicates that a 
determining element in the NN potential is its on-shell content. Indeed, we have tested this 
by taking a restricted energy range (0 - 400 MeV) for the NN phase-shifts to construct the 
inversion potential. In this case, the off-shell content of the inversion potential is different 
as reflected by the differences observed in the NA scattering observables, more significantly 
above 200 MeV. We have also considered a 0-3 GeV energy range for the genuine Paris 
phase-shifts to construct the inversion potential. In this case we observed no further differ- 
ences with respect to those observed when using the 0-1.3 GeV energy range. These results 
indicate that NN potentials which closely accounts for the NN phase-shift data over a wide 
energy range one would be able to assess the level of completeness of the optical model for 
NA scattering. 

The relatively weak sensitivity to the off-shell effects suggest that NN potentials, con- 
structed by means of the inversion scheme and following closely the NN data are meaningful 
and would provide predictions for the full-folding model very close to what would it be 
obtained from first principle NN potentials with comparable fit to the data. 



B. Effective interactions off shell 

To illustrate the degree of sensitivity of the g matrix upon alternative choices of bare 
NN potentials, we have calculated selected matrix elements relevant for the leading con- 
tribution to the optical potential. Thus, we consider diagonal g— matrix elements, i.e. 
(K,\g^ a . kF (u ; R)\k), with Q Q fixed to a single value by Q Q = \ / 2mu. In the context of 
the free t matrix (kp = 0) this kinematics implies the on-shell relative momentum occurs 
for k = \Q - Since k is in general independent of both Q Q and lu, the resulting function 
corresponds to off-shell elements of the g matrix. For the Paris and its inversion equivalent 
potential we have solved Eq. ( pf) for g in the cases /cf=0 fm _1 and kp= 1 fm -1 and for the 
state 1 Sq. This is a good example to show the differences generally observed in most of the 
states for the effective interactions. In Fig. [| we show the corresponding matrix elements, 
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both real and imaginary components, for lo= 30, 200 and 400 MeV. The solid and the dashed 
curves represent results from the Paris and the Paris inversion potentials respectively. 

Since the Paris potential and its inversion are equivalent on-shell, the solid and dashed 
curves for the t matrix (kp=0) must intercept each other at k — \Q Q in both their real 
and imaginary components. This fact is indeed the case in the figures shown here as the 
on-shell constraint is explicitly built in the inversion method; this is not necessarily the case, 
however, for the g matrix as the propagator differs from that in free space by the presence of 
Pauli blocking and the self-consistent fields. It is interesting to note the different asymptotic 
behavior (k — > oo) between the Paris and the inversion potential for both kp=0 and kp=l 
fm _1 . Whereas all curves tend to zero as k increases, the one corresponding to the real part 
of the g matrix for the Paris potential does not. This result is consistent with the explicit 
momentum dependence introduced in the parameterization of the Paris potential, feature 
which becomes dominant at high momenta. The extent to which this distinctive behavior 
is significant in the dynamics of the collision of the projectile with the nucleus needs to be 
assessed in the context of the elastic scattering observables. 

In Fig. [| we also observe that the plotted matrix elements usually differ utmost in a few 
percents around the on-shell values. This result is consistent with the differences we have 
already discussed at the level of the NA observables. 



IV. NA SCATTERING FROM NN PHASE-SHIFTS 

In this section we extend the idea of NN potentials obtained directly from NN phase- 
shifts through the quantum inversion method and its application to NA scattering. For this 
purpose, we have calculated a NN inversion potential based on the SM94 phase-shift analysis 
of Arndt and collaborators ]17|]. The main problem we face is the choice of a meaningful 



energy range where the phase-shifts are to be taken from. Indeed, in order to be consistent 
with the inversion scheme, a set of real phase-shifts are needed. This sets a limit for the 
energy range at the pion production threshold (~300 MeV). However, our studies with the 
Paris potential in sect. Ill show that phase-shifts at much higher energies are required by 
the inversion method. 

Our approach is as follows. In order to account approximately for the trend with energy 
that phase-shifts have above pion production threshold, we have neglected the imaginary 
component of the phase-shifts and retained only their real components. Thus we can con- 
struct a real NN inversion potential from a set of real phase-shifts. The range of energy 
considered for the inversion is 0-1.3 GeV. Although this approach becomes essentially qual- 
itative, we expect that our calculations will provide a guidance on how the medium and 
short range parts of the NN interaction, as determined by more realistic phase-shifts, affect 
the description of NA scattering in the intermediate energy region. 

We have calculated NN inversion potentials from the SM94 data for all NN channels with 
J < 2. In Fig. [5] we present the phase-shifts as a function of the energy for some selected 
channels. Dots correspond to the SM94 data (only the real part of the phase-shifts) and full 
curves are the results from the inversion potentials constructed from the SM94 data. We 
also plot the phase-shifts obtained from the Paris inversion potentials as reference (dashed 
curves) . These results suggest two observations. One is related to the ability of the inversion 
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method to reproduce the input data. The level of accuracy obtained is again excellent, as 
in the Paris case discussed in sect. III. The other point of physical significance is the clear 
departure of the Paris phase-shifts from, at least, the real part of the experimental ones 
above pion production threshold. It is precisely the presence of these differences which 
should affect the overall behavior of the inversion potential and, in particular, its off-shell 
components. 

The role of the NN inversion potentials based on the SM94 data has been tested in 
NA elastic scattering. We have performed in-medium (g— matrix) full-folding calculations 
for proton scattering on 40 Ca at several energies. In the actual calculations, the inversion 
potentials were used for the NN channels with J < 2. For simplicity, we have kept the 
components corresponding to the Paris potential for all the states with J > 2. In Fig. || we 
show the differential cross-section and analyzing power (A y ) results for p + 40 Ca at 40 MeV. 
The full line corresponds to the results from the inversion potential and the dashed curve to 
results from the Paris inversion one. In this and following figures we include as reference the 
calculations done with Paris inversion potentials to have a clear comparison of the differences 
that may be observed. In the scheme we have followed, these differences would come directly 
from the different underlying set of NN phase-shifts used. The results in Fig. |6] at 40 MeV 
show, however, little difference. This is consistent with the fact that both potentials are 
comparable in their agreement with the NN data below pion production threshold (see Fig. 
H). Still, some small sensitivity of the order of 5 % in the observables is observed. Since 
both potentials are constructed following exactly the same inversion procedure, we conclude 
that phase-shifts above 300 MeV still determine properties of the NN potentials which affect 
low-energy NA scattering. These small, mostly off-shell, differences stay in even at much 
smaller energies in many-body systems. 

We have pursued these calculations at higher energies. In Fig. [7] we present the results 
for p + 40 Ca at 200 and 300 MeV. The meaning of the curves is the same as in Fig. |6] and the 
data was taken from Ref. |H| at 200 MeV and from Ref. [32] at 300 MeV. At these energies 
we start noting a marked departure in the predictions of the two inversion potentials, with 
a tendency of the potential constructed from the NN phase-shift data to be relatively closer 
to the NA scattering data. This result reflects the disagreement existing between the data 
and the Paris potential phase-shifts (Fig. [|). Our findings are confirmed with calculations 
for p + 40 Ca at 400 and 500 MeV. These results are shown in Fig. [8]. Here, the data was 
taken from Ref. [3^] at 400 MeV and from Ref. |30| at 500 MeV. Certainly, applications of 



the Paris (or its inverse) potential at 400 MeV and above constitute an extrapolation of the 
model. Nevertheless, these applications serve us to illustrate both the role of the NN phase- 
shifts above pion production threshold in determining the NN potential and the ability of 
the inversion method to capture that physics. The differences given by the two potentials in 
Fig. H are remarkable. In particular, it is notable the improvement obtained in describing 
the NA scattering observables with the inversion potential constructed from the SM94 data, 
mainly for q >1 fm _1 in both da/dQ and A y . Furthermore, uncertainties associated with 
the off-shell behavior of the NN potential are smaller than the departure of the inversion 
potentials from the Paris to the SM94 data. This indicates that the improvement in the 
description of the NA data is a direct consequence of an improved account for the NN data, 
mainly above pion production threshold, a built-in feature in the inversion potential. 
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V. SUMMARY AND CONCLUSIONS 



In this paper we have addressed the problem of how properties of the underlying NN 
bare interaction determines the dynamics of a many-nucleon system. Our approach for the 
NN force is based on the quantum inversion method where a local, static, channel dependent 
potential is constructed directly from the NN phase-shifts. By departing from NN potentials 
derived from a field theoretical approach and empirically modified to maximize the fit to 
NN data below pion production threshold, we expect to shed new light on these properties 
of the NN force relevant to many-body processes. 

We have investigated effects associated with the genuine off-shell behavior of bare NN 
potentials which are equivalent from the point of view of the phase-shifts over a wide energy 
range. The generation of on-shell-equivalent potentials was based on the Gelfand-Levitan 
and Marchenko inversion method for the NN system. Applications were made with the Paris 
potential, where the inversion method was applied to its corresponding phase-shifts up to 
1.3 GeV kinetic energy in the NN laboratory system. This fairly large energy range was 
required to have the off-shell behavior of the inversion potential being determined mainly 
by the theory and not by the choice of a particular set of phase-shifts. The corresponding 
NN effective interactions were found to exhibit sizable differences off-shell, particularly for 
relative momenta above 3 fm -1 in the NN system. The investigation of these differences was 
made in the context of NA elastic scattering with the calculation of in-medium full-folding 
optical potentials for proton scattering from 40 Ca and 208 Pb and at beam energies between 
40 and 400 MeV. We found a weak sensitivity of the NA scattering observables to the off — 
shell differences observed between the NN effective interactions provided the underlying bare 
NN potentials are equivalent on-shell. These differences, at most 10% at the level of the 
NA scattering observables, are present in the whole energy range. One striking aspect is the 
ability of the inversion method to generate an NN potential which is, essentially, physically 
equivalent to its original counterpart although, by construction, their analytic properties 
and asymptotic behavior are very different. 

Based on the success of the inversion method, we have constructed inversion NN poten- 
tials based on the SM94 phase-shift analysis. These potentials have been applied to the 
calculation of full-folding optical potentials for proton elastic scattering from 40 Ca in the 
40-500 MeV energy range. We obtained results for the scattering observables which yield 
a fit to the NA scattering data at 40 MeV comparable to that obtained with the Paris po- 
tential. However, a departure from the Paris results appears above 200 MeV, with a clear 
improvement of the NA scattering data being achieved in the 400- and 500-MeV applica- 
tions. This improvement comes as a direct consequence of the closer agreement between the 
inversion potential and the NN phenomenology. An important conclusion emerging from our 
studies is that bare NN potentials which only provide a fit to NN phase-shifts below pion 
production threshold and disregard their higher energy behavior are unlikely to be realistic 
candidates to describe nucleon-nucleon dynamic in the nuclear medium. 

Despite the improvements in the description of NA scattering using inversion potentials 
from NN data, difficulties still remain in describing the NA scattering data at momentum 
transfers below 1 fm -1 . We believe the explicit treatment inelasticities of the NN interaction 
above pion production threshold and baryon excitation mechanisms need to be addressed as 
the inversion potentials from the SM94 analysis were restricted to the real part of the phase- 
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shifts only. On the other hand, in the lower energy applications we observe a systematic 
inability of the full-folding model to describe in greater detail the NA scattering data. Here, 
the presence of higher order processes due to the many-body nature of the problem need 
definitely be accounted for. Indeed, the results obtained at these energies using the Paris 
potential, its inversion and the SM94 inversion provide essentially the same description of 
the NA scattering data. Therefore, we cannot attain these limitations of the optical model 
to the uncertainties associated with the off-shell behavior of the NN interaction. 
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FIGURES 



FIG. 1. L < 2 channel phase-shifts from the Paris potential (crosses) and from corresponding 
inversion potential (full curves). 

FIG. 2. Scattering observables calculated from the Paris (solid curves) and its corresponding 
inversion potential (dashed curves) for p + 40 Ca at 40 MeV. 

FIG. 3. Scattering observables calculated from the Paris (solid curves) and its corresponding 
inversion potential (dashed curves) for p + 208 Pb at 400 MeV. 

FIG. 4. Behavior of the diagonal elements of the g matrix for the state 1 Sq as a function of the 
relative momentum and for cj=30 MeV (top), w=200 MeV (center) and cj=400 MeV (bottom). 
The figures on the left correspond to kp = ( t matrix) and those on the right to kp = 1.0 fm^ 1 
The solid curves represent results from the genuine Paris potential. The dashed curves correspond 
to results obtained from inversion potentials based on the Paris phase-shifts. 

FIG. 5. L < 2 channel real phase-shifts from SM94 data (dots), and from the inversion poten- 
tials constructed from the real part of the SM94 data (solid curves) and the Paris potential (dashed 
curves) . 

FIG. 6. Calculated and measured differential cross-section and analyzing power for p+ 40 Ca 
elastic scattering at 40 MeV. The solid and dashed curves were obtained from full-folding using the 
SM94 and Paris inversion, respectively. All curves represent in-medium full-folding calculations. 

FIG. 7. Calculated and measured differential cross-section and analyzing power for p+ 40 Pb 
elastic scattering at 200 and 300 MeV. The curve patterns follow convention of Fig. ||. 

FIG. 8. Calculated and measured differential cross-section and analyzing power for p+ 40 Ca 
elastic scattering at 400 and 500 MeV. The curve patterns follow convention of Fig. ||. 
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